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Abstract 

Assuming that the cosmological principle holds, Maartens, Ellis 
and Stoeger (MES) recently constructed a detailed scheme linking 
anisotropies in the cosmic background radiation (CMB) with anisotropies 
and inhomogeneities in the large scale structure of the universe and 
showed how to place limits on those anisotropies and inhomogeneities 
simply by using CMB quadrupole and octupole limits. First we indi- 
cate and discuss the connection between the covariant multipole mo- 
ments of the temperature anisotropy used in the MES scheme and 
the quadrupole and octupole results from COBE. Then we introduce 
those results into the MES limit equations to obtain definite quantita- 
tive limits on the complete set of cosmological measures of anisotropy 
and inhomogeneity. We find that all the anisotropy measures are less 
than 10"** in the case of those not affected by the expansion rate H, 
and less than 10~® Mpc^^ in the case of those which are. These 
results quantitatively demonstrate that the observable universe is in- 
deed close to Friedmann-Lemaitre-Robertson- Walker (FLRW) on the 
largest scales, and can be adequately modelled by an almost-FLRW 
model - that is, the anisotropies and inhomogeneities characterizing 
the observable universe on the largest scales are not too large to be 
considered perturbations to FLRW. 
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1 Introduction 

People have long considered that the almost isotropy of the cosmic mi- 
crowave background radiation (CMB) indicates that the universe on very 
large scales is isotropic and homogeneous, or nearly so. The first rigorous 
result supporting this conjecture was the theorem proved by Ehlers, Geren 
and Sachs (1968) (EGS): "If a family of freely falling observers measure self- 
gravitating background radiation to be everywhere exactly isotropic, then 
the universe is exactly Friedmann-Lemaitre- Robertson- Walker (FLRW)." 
Later, Grishchuk and Zel'dovich (1978), analyzing the FLRW perturbations, 
argued that, if CBR anisotropics are very small, all anisotropics and inhomo- 
geneities on scales larger than the horizon should also be very small. Much 
more recently Stoeger et al. (1995) (SME) proved a significant generalization 
of this theorem: "If the Einstein-Liouville equations are satisfied in an ex- 
panding universe, where there is present pressure-free matter with 4- velocity 
vector field {uau"- = — 1) such that (freely propagating) background radi- 
ation is everywhere almost-isotropic relative to u"", then spacetime is almost 
FLRW." These results obviously provide a fundamental link between ob- 
servational and theoretical cosmology - one which promises to reveal key 
aspects of the very large scale structure of the universe without relying on 
the often uncertain observational measurements of local and intermediate 
scale structures. It should be noted that both the EGS and the SME results 
depend on assuming that the cosmological principle holds - here expressed 
in terms of the almost-isotropy of the background radiation relative to every 
fundamental observer in the spacc-timc. Goodman (1995) and Maartens et 
al. (1995a), however, have recently pointed out the important fact that the 
cosmological principle itself is partially testable via the Sunyaev-Zel'dovich 
effect. 

Employing this connection between CBR anisotropics and the full range 
of cosmological isotropics and inhomogeneities, Maartens et al. (1995a,b) 
(MES) have developed a detailed scheme demonstrating how limits on the 
temperature anisotropics of the CMB imply rigorous limits on the anisotropy 
and inhomogeneity of the universe. In this paper we show how CMB 
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anisotropy data is inserted into the theoretically derived limits of this scheme 
to constrain strongly the vorticity, shear, spatial gradients, Weyl-tensor com- 
ponents, and other measures of deviation from FLRW for the observable uni- 
verse. Wc then introduce the limits COBE places on the dipole, quadrupole 
and octupole of the CBR to determine those limits quantitatively. 

In the next section we shall briefly give the key MES equations which 
constitute the CMB limits on the anisotropy and inhomogeneity of the uni- 
verse. Then, in the third section, we shall discuss the relationship of the 
usual CMB multipole results with the temperature anisotropy multipoles 
in the MES equation. Finally, in section 4, we shall present the COBE 
quadrupole and octupole observational results, and use these to arrive at 
limits on all the cosmological anisotropy and inhomogeneity measures men- 
tioned above. 

Maartens, Ellis and Stoeger (1996) recently provided such limits by in- 
troducing order-of- magnitude COBE results of £2 « €3 f« 10~^, but they did 
not discuss the relationship between temperature anisotropy multipoles used 
in their scheme and those in terms of which the COBE results are given. 
Here we fill that gap, and then use the up-to-date COBE results for the rms 
quadrupole and octupole to obtain improved limits. 

2 The MES Equations 

MES have written the covariant temperature isotropy multipoles as Tai...aL^ 
where L give the multipole number. Thus, for instance, r^. Tab, Tabc are, 
respectively, the components of the dipole, quadrupole, and octopole of the 
CMB temperature anisotropy. This form of the harmonic decomposition, 
which is given in detail in Maartens ct al. (1995a), is equivalent to that 
formulated in terms of spherical harmonics (Ellis, Matr avers and Treciokas 
1983; Ellis Treciokas and Matravers 1983). 

MES have assumed that there are observed bounds on these temperature 
anisotropy multipoles, so that there exist 0[1] constants €l such that 

I Ta,...aL \< CL- (1) 

The absolute-value brackets have been defined to be the squareroot of the 
sum of the squares of the components of a given vector or tensor (Maartens 
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et al. 1995a, p. 1526). Thus, ei gives the hmits on the dipole components, 
€2 on the quadrupole components, €3 on the octupole moments, and so on. 



Then, using the strong observational assumptions on the spatial gra- 
dients and time-derivatives of the temperature harmonics they adopted in 
Maartens et al. (1995b), they have the observational limit equations on 
various kinematic, dynamic and geometric indicators of anisotropy and in- 
homogeneity (Maartens et al. 1995a,b): 

'^"^'=4^#^<lf(12ei + ^62), (2) 



I ^ab I 5 3 

< ^ei + 362 + =63, (3) 



e 3 7 

I ^ab I 10 ^2 

< + —62, (4) 



e 3 15 



' < -He2 + (7^)[60ei + 134e2 + 663] + {-^)H[lQei + —62], (5) 



2 S7m 3 



e ' 3 



I < + ^62 + f 63) + ^{imR + 150M)i^e2, (7) 



e ' 3 3 7 45 



^afel ^ „.16 ^52 1 

< ^(^ei + —€2 + —63). (8) 



e ' 3 15 21 

Here ji and p arc the radiation and matter densities, respectively, and 
and VIm arc the ratios of the radiation energy density and the matter den- 
sities, respectively, to the critical density of the universe. H, of course, the 
Hubble parameter, and T is the CMB temperature, aah is the shear of the 
congruence of timelike geodesies in the universe, ojab its vorticity, = u?^^ its 
expansion scalar, and Eab and Hab the electric and magnetic components, 
respectively, of the Weyl tensor. These quantities, along with the spatial 
gradients of p, ji and describe the anisotropy and inhomogeneity of the 
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space-time. The Va operator expresses the covariant spatial gradient in the 
spacetime. It is defined by 

VcQa...b = hiKhl Vd Qe...f, 

where hab = dab + '^a'^b is the projection tensor in the rest spaces of the 
geodesically moving observers and Va is the covariant derivative define by 
the metric gab- Finally, it should be noted that equations (2) - (8) hold 
independently of the statistics of the underlying matter-density fluctuations 
- they do not assume that the fluctuations are Gaussian. 

If we can determine ei, €2 and €3 from CMB measurements, and have 
observational estimates of H, VLr and f^M from other observations, then we 
can determine limits on all these anisotropy and inhomogeneity indicators. 
We shall proceed to do this in Section 4. But first, in the next section, we 
need to discuss the relationship between the | Ta^...a^ | given in the equations 
and the multipole moment results determined by the COBE Differential 
Microwave Radiometers (DRM) and other CMB anisotropy detectors. 

3 CMB Multipole Anisotropy Data 

We have mentioned above that the harmonic decomposition represented by 
the I Ta^,,,aL I is equivalent to that in terms of spherical harmonics. But 
the question then is whether or not the multipole results recovered from, 
say, COBE DRM data can be simply substituted into our limit equations. 
Are the multipole results presented in the COBE papers measurements of 
the I Taj...aj;, I in our limit equations? The answer to this question is 'Yes,' 
as long as we use the real rms dipole, quadrupole, octupole moments they 
obtain, and not those associated with obtaining the power spectrum, which 
are often the focus of their reported results, determine the numerical factors 
relating these moments, defined in terms of Legendre polynomials, to the 
I Tai...aL I (see section 5 below), and as long as we realize that they are 
usually given as the squareroot of the sum of the squares of the (2L + 1) 
components of the L-pole, the rms L-pole - not as values of each separate 
component of the multipole in question. Furthermore, in the COBE data 
the multipoles are quoted for ST, instead of for which our | Ta^ . a^ | 
are the multipoles. The COBE rms values published all have units of ^K, 
therefore. To translate these values into what we need we must thus divide 
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them by T, the average CMB background temperature over the sky. 

As a relevant example, Bennett el al. (1994) give results for the square 
of the rms quadrupole amplitude Q^msj 

QLs = ^(^Q? + Ql + Ql + Ql + Ql), (9) 

where the five quadrupole components Qi arc given in terms of Galactic 
coordinates / and b (Galactic longitude and latitude respectively) by the 
expansion 



Q{l,b) = Qi(3sin^6- l)/2 + Q2sin26cosZ + Q3sin26sinZ + 

Q4 cos^ b cos 21 + cos^ b sin 21. (10) 

The 'strange' coefficients in equation (9) are due to the fact that the basis 
vectors are orthogonal but not orthonormal (see Partridge 1995); Qrms does 
not contain a factor of (2L + 1)~^ = 1/5. ^Prom four years of data the 
COBE workers give us a best fit value of (Smoot, private communication; 
Kogut et al., 1996) 

Qrms = 10.7 ± 7fiK, (11) 

with a 95% confidence limit. We can modify this result (see below) for use 
in equations (2) to (8). This quantity is to be carefully distinguished in the 
COBE results from Qmis-PS which is often referred to and which is not the 
true best fit value of the quadrupole, but rather the value of the quadrupole 
derived from a power-spectrum fit (when a power law is assumed) of the 
other higher-order multipole moments [see Bennett et al. 1994; Smoot et al. 
1992] 

The dipole moment can be neglected, according to the fairly well justified 
assumption that it is all due to our peculiar motion with respect the rest 
frame of the microwave background. In fact this is what is done in the 
COBE anisotropy analysis (Bennett et al. 1994). However, we should be 
aware that there could in principle be a small non- Doppler contribution 
to the CMB dipole (see Maartens et al. 1995a and Maartens et al. 1996 
concerning this). Thus, in our calculations below we shall set 

ei = 0. (12) 
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The COBE workers have not yet published the rms octupole from their 
data, bTit G. Smoot (private communication) has kindly informed us that 
the rms octupole results from COBE are: 

Orms = 16 ± 8fiK. (13) 

This is consistent with the results given by Wright et al. (1994) in their 
Figure 1. 

There are several important observational and data-reduction issues, and 
one theoretical issue, which we should briefly mention here, to provide the 
background against which we can understand and appreciate these COBE 
multipole results. 

The first is that the COBE DMR experiment docs not directly measure 
the dipole, quadrupole, octupole moments of the temperature anisotropy, 
but rather the two-point correlation function of the temperature anisotropy 
(Smoot et al. 1992, Padmanabham 1993, Partridge 1995) 

oo 

C(a) =< 5(n)5(m) >= ^ AT^W{efPi[cos a], (14) 

e=i 

where a is the angle between the two points, n and m are unit vectors 
denoting the two different directions, so that cos a = n-m, and the averaging 
brackets signify the average over all pairs of points on the sky with separation 
angle a. Furthermore, ^(n) = 5T{n), the temperature anisotropy in a given 
direction n on the plane of the sky. The AT^ (with £ = L) are the squares 
of the rotationally invariant rms multipole moments (thus, ATI = <3rms ~ 
see Bennett et al. 1994, for instance) 

where the aim are the coefficients of the expansion of the temperature 
anisotropy in spherical harmonics, i. e. 

S{n) = S{9, ^) = J2 (^iruYimiO, </>), (16) 

£,m 
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9 and being, of course, the angular coordinates of the direction at which 
the temperature anisotropy is being measured. In equation (14), finally, 
P£[cosa] is just the Legendre polynomial of degree £ given as a function 
of cos a, and W{i) is the window function, which describes the smoothing 
properties of the instrument's beam. For detectors measuring large-angle 
CMB anisotropies it essentially weights the multipole moments in such a 
way that the higher multipoles are smoothed over - that is, the instument 
is insensitive to anisotropies on angular scales less than a certain ^-polc, 
and W{£) describes that insensitivity and resulting transfer of power in the 
measurement from higher multipoles to lower ones. Thus, when the actual 
quadrupole or octupole is determined from the data, AT^ must be decon- 
volved from W{i). When the rms dipole, quadrupole, octupole results are 
given by the COBE researchers, this deconvolution has already been per- 
formed. This is one of the procedures which must be effected to give us the 
real rms multipole moments. 

The second major observational and data-reduction problem is that in 
the COBE measurements there is a great deal of contamination by exper- 
imental systematic errors, including Galactic emission (see Bennett et al. 
1994 and Smoot et al. 1992, and references therein) And the lower multi- 
pole moments are the most susceptible to these distortions. Furthermore, 
when the data is processed, the entire region containing the Galaxy is re- 
moved from the data set. This 'Galactic cut' destroys the orthogonality of 
the spherical harmonics and leads to further aliasing of higher order mul- 
tipole power onto the lower multipoles (dipole, quadrupole, octupole, etc.). 
Corrections for this are estimated on the basis of Monte Carlo simulations 
(Bennett et al. 1994) and included in the published values for the rms 
multipoles. There are a number of other complex issues which it has been 
necessary to resolve in arriving at these values ( see Bennett et al. 1994, 
Smoot et al. 1992, and Wright et al. 1994 for further discussion). 

Finally, there is the theoretical-observational issue of 'cosmic variance' 
Actually, cosmic variance does not affect what we are concerned with here, 
as we shall see. But it is important to realize why it does not. It may 
explain why the multipoles we measure have the values they have relative 
to theoretical models of the perturbation spectrum, but it is docs not lead 
to observational errors, which would have to be corrected for. It does af- 
fect the comparison of the measured multipole power spectrum with the 
theoretical power spectrum predicted from, say, inflationary models (Ab- 
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bott &: Wise 1984). If the primordial perturbation spectrum originated due 
to fluctuations in the inflaton field, as we think it did, then the values of 
the temperature anisotropy multipoles they induce will be random variables 
with a certain distribution, probably Gaussian with zero mean (Liddlc & 
Lyth 1993), and thus with a certain variance. Our observable universe is 
only one realization of that ensemble of universes represented by the prob- 
ability distribution. Therefore, the value for each multipole we obtain from 
our observations will give us just one point of the distribution, which, in 
general, will not reflect the ensemble averaged value. It will deviate from 
it by a certain amount, which can be theoretically estimated by the vari- 
ance (Liddle k. Lyth 1993). This variance goes as 2/{2i + 1) and so will be 
more significant for the lower multipoles (Abbott & Wise 1984; Smoot et 
al. 1992). Our concern here, however, is not to compare the observed power 
spectrum of CMB anisotropics with the theoretical spectrum of density per- 
turbations generated by an inflationary scenario. It is merely to use the best 
values of the CMB multipole moments we have available - however they are 
generated and whatever their spectrum - to set definite limits on the large 
scale anisotropy and inhomogeneity of the observable universe itself. Thus, 
cosmic variance falls outside those issues which we need to consider in ar- 
riving at those limits. 



4 Calculating the Anisotropy and Inhomogeneity 
Limits 

We are now ready to use the values for the rms dipolc, quadrupolc and 
octupole the COBE team have so far obtained to determine the anisotropy 
and inhomogeneity of the universe on very large scales. As indicated above, 
we set the dipole equal to zero - equation (12). Thus, we shall set ei = 
in our equations (2) to (8). In order to transform equations (11) and (13) 
into values of | Tab \ and | Tabc \ ■, respectively, we need to divide the results in 
equations (11) and (13) by T = 2.73^, since our multipoles are for ST/T. 
We do not need to divide them by (2£-|- 1)^/^, since our multipole quantities 
are given as the absolute values, which we defined as the squareroot of the 
sum of the squares of the components (see above). Finally, we also need to 
relate the ^Tf, the coefficients in the usual Legendre polynomial expansion, 
equation (14) to the | Tai...a( |- The numerical relationship is (Gebbie &; Ellis 
1998) 
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m 

'2^(^!) 

This gives 



r^,r^^ = (2£ + l)-y^Ar^ (17) 



I Tab P= 7.5Q2^„ (18) 

and 

I Tabc ?= 17.502^„ (19) 

We obtain, therefore, 

(ea) = 1.1(±0.8) X 10-^ (20) 

and 

(es) = 2.5(±1.3) X 10"^ (21) 



In examining equations (2) to (8), we notice that they are not expressed 
in terms of average values of the components. Let us now consider them to 
be equations for the rms averages of the indicators they represent, which we 
can do, since our e's are all positive. Now writing the Hubble parameter as 
H = 100^ km/ sec/Mpc, QA <h < 1.0, and neglecting terms containing the 
factor Oij, since this is presently so small - = 4.11^"^ x 10~^ (Roos, 
1994) - we obtain: 

^IVafA^ < 1.8/1 X 10-s Mpc-\ (22) 
A* 

(^) < 4.4 X 10-^ (23) 
(^) < 1.5 X 10-6, (24) 
/ 1 VaPl ^ < (0.02 + 0.54nj})h X 10-6 Mpc-\ (25) 



e 

) < S.Oh X 10-^ Mpc-S (26) 



e 
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( 



) < (1.5 + 480m)/j X 10"^ Mpc 



-1 



(27) 



e 



( 



Hgb 

e 



) < 1.6/1 X 10"^ Mpc 



-1 



(28) 



In equations with H we have put a factor of c back into the equations to 
give the correct order of magnitude and the correct dimensions. This factor 
is hidden in the normahzation of the 4- velocity u"", which figures imphcitly 
in the equations. 

These are the hmits we desire on the anisotropies and the inhomo- 
geneities of the universe. In the past other groups have put Umits on the 
shear and the vorticity using hmits on CBR isotropies (Bajtlik et al. 1986; 
Martinez- Gonzalez & Sanz 1995; Barrow et al. 1985). But in doing so they 
assumed exact spatial homogeneity, which yields a limit on shear which is too 
strong (Maartens et al. 1996). Our limiting scheme, as we have mentioned, 
does not assume exact homogeneity, or even that the inhomogeneities and 
anisotropics are small. The result that they are small is due to the fact that 
the CBR anisotropies are small. The primary assumption we have made is 
a weak form of the Copernican principle - that all fundamental observers 
in the relevant space-time domain measure at most the same level of CBR 
anisotropy. This, as we have pointed out, is at least partially testable - and 
fully falsifiable. A single observation demonstrating that CBR anisotropies 
relative to some other cluster of galaxies is large compared to those we ob- 
serve would banish that assumption. Finally, our approach provides limits 
on the full range of possible isotropics and inhomogeneities, including spa- 
tial gradients of the radiation and matter densities, and of the expansion 
parameter, and the electric and magnetic components of the Weyl tensor. 
It is worth pointing out that the Weyl tensor components measure those 
parts of the gravitational curvature which are not determined locally by the 
mass-energy distribution(that is, via the Einstein field equations). Instead, 
they are determined by the Bianchi identities, and so, in a sense, are due 
to the mass-energy distribution at other, more distant points (Hawking &; 
Ellis 1973; Maartens et al. 1996). Hab, the magnetic Weyl tensor, mea- 
sures the amount of gravitational radiation in the space-time. It vanishes 
in the case in which exact spatial homogeneity and isotropy are assumed. 
Eab, the electric Weyl tensor, which also vanishes in that case, measures the 
tidal, shear-inducing force of the global gravitational field. It thus manifests 
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its presence by the shear it introduces in timehke and nuU congruences of 
geodesies - worldUnes of particles. In fact, the measurement of null shear 
in bundles of light rays would be the clearest signature of the presence of 
non-zero Eab- (see Hawking &; Ellis 1973) 

We see clearly from our results, equations (19) to (25), that in every 
case the anisotropy and inhomogeneity measures for the universe on the 
largest scales are very small, despite the significant inhomogeneities which 
have been detected on intermediate scales. This provides very strong justifi- 
cation for considering the observable universe to be almost-FLRW on large 
scales - that is, the deviations from FLRW are small enough to be consid- 
ered perturbations. 

Finally, we might wonder what the relevance of these limits is to "the 
averaging problem" (see Zotov and Stoeger 1995 and references therein) in 
cosmology. Do these results demonstrate that the universe is almost-FLRW 
whatever the resolution of those issues happen to be? Unfortunately, that 
is not the case. The analysis we have implemented here assumes that the 
effective theory of gravity on cosmological scales is general relativity. The 
averaging problem calls that assumption into question. 

We are very grateful to George Smoot for providing recently derived 
COBE values for the rms quadrupole and octupole, as well as for clarify- 
ing the principal source of error and the conventions in the COBE results. 
Special thanks also to Roy Maartens and George Ellis for encouragement 
and helpful comments. Marcelo Araujo thanks the Vatican Observatory Re- 
search Group for support and hospitality during the course of this work. 
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